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(57) Abstract 

A method for continuously deteimining cardiac output (CO) in a subject comprising continuously measuring arterial blood 
pressure data at a site and converting said arterial blood pressure data to a digital signal representative of a pulse contour curve 
wayefoim, calculating heart rate, grade of slope, systolic peak pressure, the presence of reflected waves, and pressure at the dic- 
rotic notch from said digital signal; determining the area, SA, under said pulse contour curve from the start of a pulse to the dic- 
rotic notch of said pulse contour curve; for any pulse which is reflective, subtracting from said area A any portion, B, of said area 
A, representing an approximate area due to wave reflection to obtain a stroke area, SA; and calculating CO as a product of SA 
and HR. Apparatus for continuously detenmining cardiac output (CO) in a subject comprising means for cariying out each step 
of the method are also disclosed. 
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BACKGROUN D QFTHE INVENTTQN 

1. Field of the Invenrion 

This invention relates generally to a method for continuously 
determining cardiac output from a pulse contour curve based on measured 
arterial pulse data. 

2. Description of the Prior Art 

The traditional method for continuously determining cardiac output 
(CO) is the thermodilution method which requires inserting an instrument 
directly into the body of the patient. There is a safety risk with the 
thermodilution method and so others have suggested methods for 
determining CO indirectly from outside the body. Until the present 
invention, however, none of the non-invasive methods have achieved 
accuracy close enough to the accuracy of the thermodilution method. 

Warner et al, Computer-based monitoring of cardiovascular functions 
in postoperative patients. Circulation 27 (suppl 2) II 68-74, 1968, 
suggested computation of cardiac output via computer analysis of the 
central aortic pulse-wave contour. 

Warner U. S. Pat. 4,834,107 discloses a non-invasive method for 
determining heart-related parameters in patients which determines pulse 
pressure, time constant of the arterial system, systolic and diastolic 
pressure, peripheral resistance, cardiac output, and mean arterial blood 
pressure. Warner uses volume data rather than pulse data, and does not 
recognize the need to adjust for inaccuracies from reflective waves, nor, of 
course, does Warner describe a method for making the appropriate 
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adjustments. Further, the Warner method apparently must recalibrate for 
each patient, and does not provide a method for factory calibration while 
achieving very accurate cardiac output information. 

Others have suggested improvements to Warner et al's method; see, 
for example: 

English, et al. Comparison of Aortic Pulse-wave Contour Analysis and 
Thermodilution Methods of Measuring Cardiac Output during Anesthesia in 
the Dog, Anesthesiology 52:56-61.-1980; 

Wesseling et al, A simple device for the continuous measurement of 
cardiac output. Adv. Cardiovasc. Phys., 5 (part n):16-52, 1983; 

Alicandri et al. Possibility of cardiac output monitoring from the 
intra-arterial blood pressure profile, Clin, and Exper. Theory and Practice, 
A7 (2 & 3). 345-353, 1985; 

De Meersman, New Noninvasive Computerized Method for the Area 
Measurement of the Dicrotic Notch, Comput. Biol. Med., 19, No. 3:189-195, 
1989; 

Wesseling et al. Arterial haemodynamic parameters derived from 
noninvasively recorded pulsewaves, using parameter estimation. Medical 
and Biological Engineering, Nov. 1973:724-732; and 

Stem et al. Anesthesiology. 73, 3A, Sept 1990. A455. described a 
new model (pulse contour) for a peripheral arterial waveform analysis 
describing the fluid and wall motion of the arterial system. The algorithm 
is described as processing the area under the arterial waveform based on 
impedance (Zao) along the arterial tree. Stem et al appear to be using the 
Wesserling et al method, citing "A Simple Device for the Continuous 
Measurement of Cardiac Output." Adv. Cardiovasc. Phys.. vol. 5 (Part 
II):16-52. Karger. Basel 1983. Stern et al require frequent recalibrations 
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and yet they still only report achieving a coirelation of 0.65 and a 
standard deviation of ±. 0.8 liters. 

In general, the prior methods requires calibration for each patient. 
In no case is a correction for reflected pulse waves provided. 

SUMMARY OF THE INVgNTIQN 
It is an object of the present invention to provide a non-invasive 
method for accurately and continuously determining cardiac output for a 
subject without having to calibrate for said subject by use of traditional 
methods. 

It is a further object of the invention to provide an apparatus which 
provides non-invasive or minimally invasive, continuous real time cardiac 
output information from arterial pulse. 

A still further object is to provide a method for continuously 
determining cardiac output from arterial pulse measurement which 
corrects for error due to reflected pulse waves. 

These objects and others which will become apparent from the 
following disclosure are achieved by the present invention which 
comprises in one aspect an apparatus useful for continuously determining 
cardiac ouQ}ut (CO) in a subject comprising: 

means for continuously measuring arterial blood pressure data at a 
site on said subject; 

means for continuously converting said arterial blood pressure data 
to a digital signal representative of a pulse contour curve waveform; 

means for continuously calculating heart rate (HR), grade of slope, 
systolic peak pressure, the presence o: reflected waves, and pressure at 
the dicrotic notch from said digital signal; 
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means for determining the area, SA, under said pulse contour curve 
from the start of a pulse to the dicrotic notch of said pulse contour curve; 

for any pulse which is reflective, means for subtracting from said 
area A any portion, B, of said area A, representing an approximate area 
due to wave reflection to obtain a stroke area, SA; and 

means for calculating CO as a product of SA and HR, 

In another aspect, the invention is directed to a method for 
continuously determining cardiac output (CO) in a . subject comprising: 

continuously measuring arterial blood pressure data at a site on said 
subject; 

continuously converting said arterial blood pressure data to a digital 
signal representative of a pulse contour curve waveform; 

continuously calculating heart rate (HR), grade of slope, systolic peak 
pressure, the presence of reflected waves, and pressure at the dicrotic 
notch from said digital signal; 

determining the area, SA, under said pulse contour curve from the 
Stan of a pulse to the dicrotic notch of said pulse contour curve; 

for any pulse which is reflective, subtracting from said area A any 
portion, B, of said area A, representing an approximate area due to wave 
reflection to obtain a stroke area, SA; and 

calculating CO as a product of SA and HR. 

BRIEF DRSrR TPnON QFTHRDRAWINfiS 

The drawings illustrate but a few of several possible embodiments of 
the invention. 

Other objects and many cf the attendant advantages of this 
invention will be readily appreciated as the same becomes better 
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understood by reference to the following detailed description when 
considered in connection with the accompanying drawings wherein: 
FIG. 1 is normal pulse curve; 

FIG.2 is a reflective curve showing determination of areas to be 
subtracted in determining SA; 
FIG. 3 is a damped curve; 
FIG. 4 is a special reflective curve; 

FIG. 5 is another special reflective curve showing Area A; 

HG. 6 is the same reflexive curve as FIGS. 5 and 6, showing Area B; 

and 

FIG. 7 is a perspective view of a hand of a subject and a front 
elevational view of the overall apparatus of the invention. 



DETAILED DESCRIPTION OF THE INVENTION 
AND THE PREFERRED EKfBODIMENTS 

Referring now in greater detail to the various figures of the drawings 
wherein like reference numerals refer to like parts, in FIG. 7 an apparatus 
for continuously determining cardiac output (CO) in a subject 12 is 
illustrated wherein a Finapres ® device for continuously measuring 
variation in arterial blood pressure data at a finger having a finger collar 
13, wrist mounted transmission module 14, wire 15, and computer which 
outputs arterial blood pressure data 16; an analogue to digital converter 17 
for continuously converting said arterial blood pressure data to a digital 
signal representative of a pulse contour curve waveform 18. A computer 
19 programmed with an algorithm, described in detail below, functions to 
continuously calculate cardiac output (CO) from said digital signal with 
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accuracy approaching or matching that of the standard thermodilution 
method. The apparatus also preferably includes an amplifier 20, a 
preamplifier 10, and a tape storage device for digital information 11. 

As an alternative embodiment to the non-invasive Finapres device, a 
minimally invasive arterial line (not shown) for continuously measuring 
blood pressure at any peripheral site can be used. 

It is well known that CO is the product of stroke volume (SV) and 
heart rate (HR). SV is the volume the heart puts out at each pump^ Our 
method comprises assuming stroke area (SA) under the pulse pressure 
curve (PPC) between the start of a pulse (upstroke) to the dicrotic notch is 
proportional to SV. We discovered that while SA is roughly proportional to 
SV, certain very critical corrections are necessary and we have discovered 
a method for making such corrections on a continuous basis, and that 
accurate CO information can be obtained with new patients without ever 
needing calibration. This information includes the absolute CO and CO 
trends during an operation, in the intensive care unit, and for outpatients 
who do not need surgical procedures. 

Each patient has a different pulse wave, and within a single patient, 
there are pulse to pulse variations, not only in frequency, corresponding to 
HR, and amplitude, corresponding to pulse-pressure, but also variations in 
shape due to signal distortions and presence or absence of reflective waves 
and damped waves due to changes in the vascular impedance (VI). These 
variations are different at different sites of measurement, and are 
dependant on the patient's vasculature condition. We have discovered that 
CO can be determined more accurately than previous methods by 
continuously calculating beat-to-beat the VL Variations in true CO are due 
to variations in ventilation, changing hemodynamic conditions and 
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medications administered, and a continuous, correct reading is necessary 
during surgical procedures to assist the physician in continuous evaluation 

of these conditions. 

Under our method, an algorithm is provided which computes (HR) 
(also called "pulse period"), grade of slope, systolic peak pressure, diastolic 
pressure, pulse pressure, mean pressure, diastolic area, the presence of 
reflected waves, and pressure at the dicrotic notch from said digital signal, 
- and SA. The preferred method is using a decaying exponential expression 
of the form Y = A + B*exp(-K*pp). Different constants (A, B, K,) are used at 
different Mean Pressures. An expression for VI can be obtained by fitting 
a curve to data or generating a theoretical expression. VI is calculated 
beat-to-beat in the preferred embodiment because VI is never constant. 

More specifically, the algorithm defines a pulse as that portion of a 
waveform beginning with a sharply rising segment (upstroke), 37 in FIG. 1, 
which is recognized by the algorithm by its slope exceeding a preset 
critical value and said value is maintained over a preset period of time. 
The end 26 or 32 of one pulse denotes the beginning of the next. In order 
to distinguish tiie beginning of a pulse from a prominent dicrotic notch, tfie 
algorithm determines whether tiie pulse pressure between the relative 
minimum and relative maximum exceeds a predefined critical value, 
preferably 50% of previous pulse pressure. If it does exceed said value, 
tiie algorithm identifies the point as the beginning/end of a pulse. The 
metiiod distinguishes tiie beginning of a pulse from the end of a prominent 
dicrotic notch by comparing the slopes and identifying the one witii tiie 
greater slope to be tiie upstroke of the pulse. Furtiiermore, larger pulse 
pressures have steeper upstrokes and the program also accounts for that 
in distinguishing upstroke from part of tiie dicrotic notch. 
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To detennine dicrotic notch, the algorithm uses pulse period as an 
index into a table of start/end offsets which delimit the boundaries of the 
dicrotic notch and within these boundaries determines where the 
instantaneous slope exceeds a critical slope value, then assigns that point 
22 as the notch. In cases of a damped curve wherein the curve of the 
notch does not exceed the critical slope value, the end of the offset is 
assigned as the location of the notch. 

The program handles three different types of pulse wave curves 
differently: reflective (HG. 2). normal (HG. 1), and damped (HG. 3); and 
thereby corrects the SA to greatly improve the accuracy of the CO 
calculation. 

A preferred method of determining whether a reflective wave exists 
is to first find the first apparent systolic peak 21 (FIGS. 1 and 2) where 
grade of slope is 0, determine the location of the dicrotic notch 22 by 
conventional steps, divide the distance along the X-axis between the first 
apparent systolic peak 21 and the dicrotic notch 22 in half to detennine 
search limit 23. and to define a search distance 24 between the first 
apparent systolic peak 21 and the search limit 23, and then search the 
search distance 24 for a poim 25 where the slope exceeds a critical value 
and label said point 25 (HG. 2) as the true systolic peak and the original, 
apparent systolic peak 21 as a false systolic peak. 

When a reflective wave is determined (HG. 2), a line 29 is calculated 
from point 21 with a slope equal to the most negative slope between the 
time values. 38 and 31. of points 21 and 22. respectively, and a point on 
said line 31. The area 27 above line 29 is assigned as "Area B" and is 
subtracted from the area under the curve between the time of the 
beginning of the pulse 26 and the time value of the dicrotic notch 31 (Area 
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A) for purposes of calculating SA. In curves such as illustrated in FIG. 4 
where line 29 intersects a line 34 representing the diastolic pressure 
before it intersects line 31, an area 33 is defined as the area under the 
pulse pressure curve bounded by line 35 from apparent systolic peak 21 
and point 36 where Hnes 34 and 31 intersect (FIG 6 shows this area as 
Area B), and is subtracted in determining corrected SA. 

For any pulse which is reflective, the computer algoritfim subtracts 
from or adds to_ said area A any portion, B, of said area A, representing an 
approximate area due to wave reflection to obtain a stroke area, SA; and 
the algorithm calculates CO as a product of SA and HR. In the preferred 
embodiment, the algorithm corrects the resultant CO by multiplying using 
an admittance factor, Y. which is a function of pulse pressure and mean 
pressure. A relationship between Y and CO/SA-HR is deteimined as a 
function by comparing the thermodilution derived CO's for each patient in 
the studypopulation over a range of pulse pressures, mean pressures, 
heart rates, and vascular impedances. Very accurate, reproducible CO's 
and CO trends can be obtained using a dynamic admittance expression, 
exceeding the accuracy and reproduceability of other pulse contour 
methods. The Y function is preferably programmed into the algorithm so 
that corrected CO information is continuously available for new subjects 
without the need to calibrate for any subject with thermodilution methods 
or any other invasive method. 

The contribution reflective waves make to the stroke area changes 
depending on the location and position the data is recorded from. For 
example, the data from the radial artery of the wrist will differ from the 
data from the femoral artery of the groin. 



wo 92/11804 PCT/US92/00408 



When the pulse pressure wave is normal (FIG. 1) or damped (FIG. 3) 
no subtractions or additions are necessary. 

As a preferred embodiment, an algorithm is presented below in Turbo C 
language showing the preferred constants, lookup tables, and the like. 

double gen_grat2line(struct curve_data curvestuff, struct notch.data notchstuff 
int the^gap, int *reverb, int *newpeak) 

{ 

int xl,x2,yl,y2; 
int minslope = 0; 
int testpoint; 

int i, data^point, extrapolation; 
double m, b, new_sum = 0; 
int firstval, temp; 

if (the^ap >= 60) 

'''newpeak = peak_corrected(curvestuff,notchstuff); 

xl = (curvestuff). top; 
yl = buff.ptr[xl]; 

for (i = xl; i <= notchstuff. start-5; i++) 
{ y2 = buff.ptr[i+l]; 

m = y2 - yl; 

yl = y2; 

if (m < minslope) 
minslope = m; 

) 

xl "= curvestuff. top; 
yl = buff_ptr[xl]; 

x2 s (notchstuff).start; 
b a yl - (in*xl); 

testpoint = (in*(notchstuff.start)) + b; 
if (testpoint < buff_ptr[curvestuff.start]) 
{x2 » (notchstuff.start); 
y2 = bu£f_ptr[ciirvestaff.staTt]: 
m = (y2 - yl)/(x2 - xl); 
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b = yl - (in*xl); 

} 

firstval = bufLptr[curvestuff.start]; 

for (i = curvestuff.start; i <= xl; ++i) 

new^sum = new_sum + buff_ptr[i] - firstval; 
♦reverb = 0; 

for(i = xl+l;i<-x2;++i) 
{ 

data_point = biifLptr[i]; 
if (data_point > *newpeak) 

♦revCTb = *reverb + (data_point - *newpeak); 
data_poiiit - *newpeak; 

1 

extn^olation = (m*i) + b; 

if (extrapolation > data_point) 
{temp = data^point - firstval; 
if (temp < 0) 

temp = 0; 
new_simi = new_sum + temp; 

) 

else if (extrapolation > firstval) 

{new_sum = new_sum + extrapolation - firstval; 

) 

} 

return(new_sum); 

} 

char is reflection(struct ciirve_data curvestuff) 

{ " 
int start s curvestuff-top, 

end = start+10, 

i = start, 

slope; 

char done = FALSE; 

wWle (0 < end) AND (!done)) 

slope = buff_ptr[i+2] - buff_ptr[i3; 
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if (slope >= 0) 
done = TRUE; 



) 

return(done); 

) 

int peak.corTected(struct curve^data curvestuff, struct notch.data notchstuff) 
{ 

#defin6 spike -60 
int i, 

start = curvestuff.top, 
double temp; 

int slope,yl,y2; — - 
int threshold = buff_ptr[stait]; 
char done » FALSE; 

temp = notchstuff .start-5; 
temp = (start + temp)/2; 
end = temp; 
i = start+S; 
while (Idone) 

{ 

yl = buff_ptr[i]; 
y2 = buff_ptr[i+3]; 
slope = y2 - yl; 
if (slope > spike) 

{ 

done = TRUE; 
threshold = yl; 

) 

if (++i >=end) 
done = TRUE; 

} 

if (threshold < buff_ptr[start]) 
{threshold = buff_ptr[i]; 

) 

retum(threshold); 

) 

double get_startoffset(double period) 
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{ 

double time^offset; 
double temp; 

if (period > 1.4) 

time^offset = (double) timeconst*(0.18); 
else if (period > 0.7) 

time.offset = (double) timeconst*(0.16); 
else if (period > 0.65) 

time^offset = (double) tiineconst*(0.14); 
else if (period > 0.6) 

time.offset = (double) timeconst*(0.13); 
else if Opcriod > 0.55) " 

time.offset = (double) timeconst*(0.12); 
else if (period > 0.5) 

time_offset = (double) timeconst*(0.11); 
else if (period > 0.4) 

time^offset = (double) timeconst*(0.08); 
else if (period <= 0.015) 

time_offset = (double) iimeconst*(0.015); 
else 
{ 

temp = (period - 0.015)*(period - 0.015); 
temp = temp*temp; 
temp'= temp*temp; 

time.offset = (double) timeconst*( 0.02 + (2.5*temp) ); 

} 

retum(time_offset); 

) 

double get_Qoffset(double period) 
{ 

double Qoffset; 

if (period > 1.2) 

Qoffset = 0.45; 
else if (period > 1) 

Qoffset = 0.5; 
else if (period > 0.8) 

Qoffset = 0.58; 
else if (period > 0.7) 

Qoffset = 0.68; 
else if (period > 0.5) 
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Qoffset = 0.75; 
else 
Qoffset = 0.95; 

retum(Qoffset); 

) 

void scan_notch(struct curve^data ciirvestuff,struct notch^data *notchsiuff) 
{ 

int yl,y2,slope; 
int upcount = 0, 
downcount = 0; 
int notch_loc, 

notch_x; 
char done = FALSE; 
char done2 - FALSE; 

int critval = -5, 
crittimes = 1; 
int _x,tt = 0; 
int i; 

unsigned int start,end; 

double twave_offset = 0; 
double pulse^period, Qoffset; 

pulse_period = (double) (curvestuff.endx - curvestuff.start)/timeconst; 
twave_offset = get_startoffset(pulse_period); 
Qoffset = get_Qoffset(pulse_period); 

start = curvestuff.top + twave_offset; 
end = (int) (curvestuff.start 

(curvestuff.endx-Cttrvestttff.start)*Qoffset); 

_x = curvestuff.top_x + twave.offset; 

yl = buffjptr[startl; 
i s start; 

while ((i <= end) AND (Idone)) 
{ 

_x++; 
i-H-' 

if bi >= 640) 
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_x = 0; 

y2 = buff j)tr[i]; 
slope = y2 - yl; 
yl = y2; 

if (slope >= critval) 

{ 

upcount++; . 
if (upcount >= 2) 
downcount = 0; 

) 

else if (slope < critval) 

{ 

downcount++; 
if (downcount >= 5) 
upcount =5 0; 

}; 

if (upcount >= crittimes) 

( 

notchjoc = i - (crittiraes+1); 
notoh_x = _x - (crittimes+1); 
done = TRUE; 

) 

) 

if (done ~ FALSE) 
{ 

notchjoc = i - (crittimes+1); 
notch_x = _x - (crittimes+1); 
done = TRUE; 

1 

if (done) 
( 

(*notchstuff).start = notchjoc; 
(*notchstu£f).x = notch_x; 

} 

else 

(*notchstuff).start = 0; 

if (done) 
{ 
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while ((!done2) AND (i <= end)) 

{ 

slope = buff_ptr[i+2] - buff_ptr[i]; 
if (slope > 0) 
done2 - TRUE; 

tt++; 

) 

if (donel) 

{ 

notchjoc = i; 
notch_x = _x + tt; 

) 

(*notchstuff).staTt = notchjoc; 

(*notchstuff).x = notch_x; — 

} 

) 

void getarea(struct curve_data *curvestuff) 

{ 

int i; 

for (i = (*curvestuf0.start; i <= (*curvestuff).endx; 1++) 

(*curvestuff).suinl = (*curvestuf0.suml + buff_ptr(i] - curvemin; 



double getnotch_area(struct curve_data curvestuff^struct notch.data notchstf) 

double tempsum = 0; 
double .temp; 
int baselength,i; 
int firstval,yval, temp; 
int xx; 

XX = curvestu£f.xl; 

firstval = buff.ptr[curvestuff.start]; 

baselength = notciistf.staTt - curvestuff.stait + 1; 
for (i = curvestuff.start; i <= notchstf.start; ++i) 

{ 

if (XX >= 640) 
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XX = 0; 

yval = buff_ptr[i]; 
if (yval < firstval) 
yval = firstval; 

tempsum = (double) (tempsum + yval); 
XX++; 

^_temp = (double) baselength*firstval; 

tempsum - (double) (tempsum - ^temp); 
xetum(tempsum): 

} 

double admittance(int mean_piessure. int pulse.pressure. int notch^ap) 
( 

double Y; 

if (mBan_pressure >= 100) 
{ Y = 1.94 + (18.01*exp(-0.14*pulse_pressure)); 

if (notclugap <= 0) 
Y4=Y*(0.5); 
else if (notch_gap <= 3) 

Y +«= Y*(03); 

if ((pulse_pressute >= 75) AND (pulsej)ressure <= 79)) 
{ if (notch_gap < 10) 

notch_gap = 10; 
else if (notch_gap > 40) 

notch_gap = 40; 

Y = Y + (notch_gap - 25)*(0.0369); 

^ else if ((pulse4)iessuie >= 50) AND (pulsejjressure <= 59)) 
{ if (notch_gap < -10) 
notch_gap = -10; 
else if (notch_gap > 25) 
notch_gap = 25; 

Y = Y - (notch_jsq) - 25)*(0.0355); 

^ else if ((pulse_pressnre > 140) OR (pulsej)ressure < 40)) 

if (alarai) printf("%c",7); 
} /* end MP >= 100 */ 
else if (mean_pressure >= 95) 

{ /*Y = 2.73 + (3.25)*exp(-0.035*pulse_pressure);*/ 

Y = 0.04 + (5)*exp(-0.008*pulse_pressure); 
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if (notch_gap <= 0) 
Y4«Y*(0^); 
else if (notch_gap <= 3) 

Y Y*(0 3); 

if ((pulsejpressure > 90) OR (pulse_pressure < 45)) 
if (alarm) printf("%c",7); 
) /* end MP >= 95 */ 
else if (mean_pressure >= 90) 
{ Y = 1.4 + (16.8)*exp(-0.038*pulse_j»ressure); 

if (notch_gap <= 0) 
Y+=:Y*(0^); 
else if (notclugap <= 3) 

Y += Y*(0 3)*» 

if ((pulse jressure > 120) OR (pulse_pressure < 40)) 
if (alam) printf("%c",7); 
) /* end MP >- 90 */ 
else if (mean_pressure >= 85) 
{ Y = 0.12 + (4.3)*exp(-0.007*pulse_pressure); 

if (aotch_gap <= 0) 
Y4«Y*(0^); 
else if (notch_gap <= 3) 

Y += Y*(0 3); 

if (0>ulsejpressure > 90) OR (pulse_pressure < 40)) 
if (alann) printf("%c",7); 
) /* end MP >= 85 */ 
else if (mean_pressure >= 80) 
{ Y = 0.26 + (4.6)*exp(-0.009*pulse_pressure); 

if (notch_gap <= 0) 
Y+=Y*(05); 
else if (notch_gap <= 3) 

Y += Y*(03); 

if ((pulse^pressure > 80) OR (pulse_pressure < 35)) 
if (alarm) printf("%c".7); 
) /* end MP >= 80 */ 
else if (mean_pressure >= 75) 
{ Y = 2.50 + (10.97)*exp(-0.068*pulse_pressure); 

if (notch_gap <= 0) 
Y+=Y*(0.5); 
else if (notch_gap <= 3) 

Y+=Y*(03); ,..v 
if ((pulse_pressure >= 60) AND (pulse_pressure <= 64)) 

{ if (notch_gap > 25) 
notch^ap = 25; 
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else if (notch_gap < 5) 
notch_gap = 5; 

Y = Y - (notoh_gap - 20)*(0.1197); 

^ else if ((pulse_pressure >= 70) AND (pulse^pressure <= 74)) 
{ if (notch_gap < 10) 
notch_gap = 10; 
else if (notch_gap > 20) 
notch_gap = 20; 

Y = Y - (notch_gap - 20)*(0.077); 

^ else if ((piilse_pressure >« 50) AND (pulse_pressure <= 54)) 
{"if (notcluf^ < 5) 
notch_gap = 5; 
else if (notch_gap > 25) 
notch_gap = 25; 

Y = Y - (notch_gap - 20)*(0.0752); 

^ else if ((pulse_pressure > 85) OR (pulse_pressure < 30)) 

if (alarm) printf("%c",7); 
) /* end MP >= 75 */ 
else if (mean_pressure >= 70) 
J Y = 1.31 + (16.03)*exp(-0.05*pulse_pressure); 
if ((pulse_pressure 60) AND (pulsejressure 64)) 
{ if (notch_gap > 20) 
notch_gap = 20; 
else if (notch_gap < 0) 
notch_gap = 0; 

Y = Y - (notch_gap - 20)*(0.136); 

^ else if ((pulse^pressure >= 55) AND (pulse_pressure <= 59)) 
{ if (notch_gap > 20) 

notcb_gap = 20; 
else if (notch_gap < 5) 
nolch_gap = 5; 

Y = Y - (notch_gap - 20)*(0.1257); 

^ else if ((pulse_prcssore >= 80) AND (pulse_pressure <= 84)) 
{ if (notch_gap < 25) 
notch_gap = 25; 
else if (notch_gap > 30) 
notch^gap = 30; 

Y = Y + (notch_gap - 20)*(0.1257); 
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else if ((pulse_pressure > 85) OR (pulse_pressure < 30)) 
if (alam) printf("%c",7); 
} /* end MP >= 70 */ 
else if (mean_pressure >= 65) 
{ Y = 1.7 + (5.69)*exp(-0.028*pulse_pressure); 
if ((pulse^pressure >= 50) AND (pulse_pressure <= 54)) 

( if (notclusaP < 5) 
notch^ap = 5; 
else if (notch_gap > 30) 
notch_gap = 30; 

Y = Y - (notch_£ap - 20)*(0.0743); 

^ else if ((pulse^pressure >= 40) AND (pulse_pressure <= 44)) 
{ if (notch_gap < 4) 

notch_gap = 4; 
else if (notch_gap > 15) 
notch_gap = 15; 

Y = Y - (notch_gap - 20)*(0.0821); 

^ else if ((pulse_pressure > 70) OR (pulse_pressure < 20)) 
. if (alam) printf("%c",7); 

) /* end MP >= 65 */ 
else if (meaiupressare >= 60) 
{ Y is 232 + (20.28)*exp(-0.066*pulse_pressure); 
if ((pulse_pressure >= 70) AND (pulse_pressuie <= 79)) 
{ if (notch_gap < 15) 
notch_gap - 15; 
else if (notch_gap > 30) 
notch_gap = 30; 

Y = Y + (notclugap - 20)*(0.0456); 

^ else if ((pulse_pressure >= 45) AND (pulse_pressure <= 49)) 
{ if (notcfa^ap < 5) 
notch_gap = 5; 
else if (notch_gap > 25) 
notch^ap = 25; 

Y = Y - (notch_gap - 20)*(0.0784); 

^ else if ((pulse^pressure >= 65) AND (pulsej)ressure <= 69)) 
{ if (notch_gap < 14) 
notch_gap = 14; 
else if (notch_gap > 25) 
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notch_jap = 25; 
Y = Y + (notch_gap - 20)*(0.0519); 

^ else if ((pulse_pressure > 80) OR (pulsej>ressure < 30)) 

if (alarm) printf("%c",7); 

) /* end MP >= 60 */ 

else V 
{ Y = 2.64 + (6.16)*exp(-0.056*pulse_pressure); 

if ((pulsejressure > 70) OR (pulse_pressure < 15)) 

if (alarm) printf("%c",7); 
) /* end MP < 60 */ 
retum(Y); 

) ~" 

void do_calcs(struct curve.data curvestuff,struct notch.data notchstf. 
double .adjust) 

{ 

double notcharea; 

double totalarea = curvestuff.suml; 

int timeinterval = curvestuff.endx - curvestuff.start; 

int frequency = (60*timeconst/timeinterval); 
double mean_pressure; 
double sttokevolume^troke; 
double notchgap, 

top_notchgap; 
double temp,Y_factor; 
int top_mark, gap, newpeak; 
int newgap, reverbs; 
double CO,C02. resistance; 

int decaytemp; 
int sign2; 

char reflects = FALSE; 
double decayslope; 
double Qpressure, avgQ; 
double gapdouble; 
int mp, pp. ng; 

top_mark = curvestuff.top; 

syU^SS^eH^^u^^^^ - curvemin)/(curvemax - curvemin); 
top_mark = curvestuff.start; 
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dia = buff_ptr[top_markl; 

dia = (double) piessure*(dia - curvemin)/(curvemax - curvemin); 
gap = sys - dia; 

notchgap = buff_ptr[notchstf. start]; 
top^notchgap = buff_ptr[curvestuff.top]; 

notcheap = pressure* 

(notchgap - buff_ptrlcurvestuff.startl)/(curvemax-curveinm). 

top^notchgapj Fcssure^ buff.ptr[notchstf.start))/(curvemax-curvemin); 
notcharea = getnotch_area(curvestuff7notchstf); 

mean^pressure = totalarea/(timeinterval*_adjust); 

if (newpeak == buff_ptr[curvcstuff.top)) 
newgap = gap; 
else 

^ temp = (double) pressure*(newpeak-cuTvemin)/(curvemax-curvemin); 
temp = temp - dia; 

newgap = temp + ((0.5)*(gap-temp)); 
avgNG = avgNG + notchgap; 
reflects = is_reflection(curvestuff); 
if ( ((notchgap <= 13) OR (gap < 60)) AND ('.reflects)) 

strokevolume = notcharea; 
dse if (sys < 140) 

^ strokevolume = (strokevolume + reverbs); 
) 

else 

^ strokevolume = (strokevolume); 
gapdooble = gap; 

newgap = temp + ((0.75)*(gapdouble-temp)); 
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) 

mp = (int) mean^pressure; 
ng = (int) notchgap: 
pp s gap; 

YJactcM- = admittanceCmp, pp, ng); 

stroke = strokevolume; ^ 

stroke = stroke*Y_factor/Cadjust*nineconst); 

CO « BSA*stroke*frequency; 

} /* end procedure do.calculations */" 

Various alternatives, modifications, and improvements should 
become readily apparent to those skilled in the art reading the above 
specification without departing from the spirit and scope of the mvenuon 
as set forth in the following claims. 
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1. A method for continuously determining cardiac output (CO) in a 
subject comprising: 

continuously measuring arterial blood pressure data at a site on said 
subject; 

continuously converting said arterial blood pressure data to a digital 
signal representative of a pulse contour curve waveform; 

continuously calculating heart rate (HR), grade of slope, systolic peak 
pressure, the presence of reflected waves, and pressure at the dicrotic 
notch from said digital signal; 

determining the area, SA, under said pulse contour curve from the 
start of a pulse to the dicrotic notch of said pulse contour curve; 

for any pulse which is reflective, subtracting from said area A any 
portion, B, of said area A, representing an approximate area due to wave 
reflection to obtain a stroke area, SA; and 

calculating CO as a product of SA and HR. 

2. Method in accordance with claim 1 further comprising correcting 
the resultant CO by multiplying using an admittance factor, Y, which is a . 
function of pulse pressure, said function having been previously 
determined statistically over a patient population by comparison of the 
product of SA and HR for a particular patient with thermodilution cardiac 
output data over a range of pulse pressures for said patient. 

3. Method in accordance with claim 2 wherein said Y is a standard 
curve which is used for every subject. 



wo 92/11804 



26 



PCr/US92/00408 



4. Method in accordance with claim 1 wherein said admittance 
factor, Y, had been determined for a particular site of measurement of 
arterial blood pressure in each of the patients in said patient population. 

5. Method in accordance with claim 1 wherein said portion B is 
determined by a method which comprises continuously calculating for each 
pulse the portion of the area under the pulse contour curve from where 
grade of slope equals zero and then increases by a set-amount, 
representing the beginning of a reflected wave* to a point on said curve 
where said grade of slope remms to zero, said point where said grade of 
slope returns to zero representing the dicrotic notch. 

6. Method in accordance with claim 5 wherein said set amount is 
greater than -200 mm Hg-sec-1 and said point on said curve where said 
grade of slope returns to zero is greater than -60 mm Hg-sec-1. 

7. Method in accordance with claim 5 wherein said set amount is 
determined by a formula based on the pulse pressure of the previous 
pulse. 

8. Method in accordance with claim 5 whoein a corrected value 
representing the true systolic peak is determined for each pulse which 
contains a reflected wave and a resultant elevated, false systolic peak, and 
an extrapolation line from said true systolic peak and the dicrotic notch is 
calculated, said extrapolation line marking a lower boundary of said area, 
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B, said area B having an upper boundary marked by said pulse contour 
curve including said elevated, false systolic peak. 

9. Method in accordance with claim 8 wherein an additional area, C, 
defined as the area above the true systolic peak, is determined and 
subtracted from said area A* 



10. Method in accordance with claim 8 wherein said area B is bound 
on one side by the false systolic peak and on the opposite side by a point 
on said pulse contour curve for said pulse where the time ordinate equals 
that of the dicrotic notch or the pressure is the diastolic pressure. 

11. Method in accordance with claim 1 wherein said site is an index 
finger using a device which measures variation in blood pressure at said 
site. 

12. Method in accordance with claim 1 wherein said site is a 
radialanery and is measured by means of a catheter. 

13. Method in accordance with claim 1 wherein said Y is calculated 
according to a formula 

A + B«e-(K(mp)«pp) 

where A is a function of mean pressure, B is a function of mean pressure, K 
is a function of mean pressure, and pp represents pulse pressure. 
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14. Method in accordance with claim 1 further comprising 
continuously calculating diastolic pressure, pulse pressure, mean pressure, 
and diastolic area. 

15. Apparatus for continuously determining cardiac output (CO) in a 

subject comprising: 

means for continuously measuring arterial blood pressure data at a 

site on said subject; 

'means for continuously conv«ting said arterial blood pressure data 
to a digital signal representative of a pulse contour curve waveform; 

means for continuously calculating heart rate (HR). grade of slope. 
systoUc peak pressure, the presence of reflected waves, and pressure at 
the dicrotic notch firom said digital signal; 

means for determining the area. SA. under said pulse contour curve 
from the start of a pulse to the dicrotic notch of said pulse contour curve; 

for' any pulse which is reflective, means for subtracting from said 
area A any portion. B. of said area A. representing an approximate area 
due to wave reflection to obtain a stroke area, SA; and 

means for calculating CO as a product of SA and HR. 

16. Apparatus according to claim 13 further comprising means for 
correcting tiie resultant CO by multiplying using an admittance factor. Y. 
which is a function 

of pulse pressure, said function having been previously determined 
statistically over a patient poputetion by comparison of die product of SA 
and HR for a particular patiem with thermodilution cardiac output data 
over a range of pulse pressures for said patient. 
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17. Apparatus according to claim 13 which comprises means for 
continuously calculating for each pulse the portion of the area under the 
pulse contour curve from where grade of slope equals zero and then 
increases by a set amount, representing the beginning of a reflected wave, 
to a point on said curve where said grade of slope exceeds -60 mm Hg-sec- 
1, said point where said grade of slope exceeds -60 mm Hg-sec-l 
-representing the dicrotic notch. 

18. Apparatus according to claim 13 wherein said means for 
continuously measuring arterial blood pressure data at a site on said 
subject comprises means for continuously measuring variation in pressure 
at an finger or toe site and means for continuously converting the resultant 
analogue data to digital data representing pulse pressure. 

19. Apparatus according to claim 13 wherein said means for 
continuously measuring arterial blood pressure data at a site on said 
subject comprises catheter means for continuously measuring pulse 
pressure data at a radialartery. 
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